Non-thermal X-rays from the Ophiuchus galaxy cluster and dark matter annihilation 
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We investigate a scenario where the recently discovered non-thermal hard X-ray emission from the 
Ophiuchus cluster originates from inverse Compton scattering of energetic electrons and positrons 
produced in weakly interacting dark matter pair annihilations. We show that this scenario can ac- 
count for both the X-ray and the radio emission, provided the average magnetic field is of the order 
of 0.1 fj,G. We demonstrate that GLAST will conclusively test the dark matter annihilation hypoth- 
esis. Depending on the particle dark matter model, GLAST might even detect the monochromatic 
line produced by dark matter pair annihilation into two photons. 

PACS numbers: 95.35.+d„ 98.80.Cq, 95.85.Nv, 95.85.Pw 



Clusters of galaxies are the largest bound dark matter 
(DM) structures in the universe. As such, they are nat- 
ural targets for the search for observational signatures 
of particle DM If DM is in the form of weakly in- 
teracting massive particles (WIMPs) 0, Q , DM pair an- 
nihilations generically produce 7-rays as well as a non- 
thermal energetic electron-positron (e^) population. The 
latter, in turn, is expected to yield secondary emissions at 
soft 7-ray, X-ray and radio frequencies via inverse Comp- 
ton scattering, bremsstrahlung and synchrotron radia- 
tion, opening up the possibility of a multi-wavelength 
approach to particle DM detection [3, Q] . A generic fea- 
ture of the broad-band DM annihilation spectrum is a 
significant hard X-ray component [l|, [B[ . 

Interestingly, the discovery of a non-thermal hard X- 
ray emission from the Ophiuchus cluster, detected with 
relatively robust statistical significance in a 3 Ms obser- 
vation with the IBIS/ISGRI and JEM-X instruments on 
board INTEGRAL, was recently reported in Ref. Q . The 
Ophiuchus cluster is a nearby (z ~ 0.028 Q) rich cluster 
with a high temperature plasma (kT ~ 10 keV), featur- 
ing the second brightest emission in the 2-10 keV band. 
In addition to X-ray observations, the steep-spectrum ra- 
dio source MSH 17-203 [8J] was associated to the Ophi- 
uchus cluster 0], indicating the presence of relativistic 
electrons. However, the cluster was not detected at 7- 
ray frequencies by EGRET [9j . 

Several clusters are known to host extended radio 
emissions [loj ]. suggesting the existence of energetic 
non-thermal electrons that radiate at radio frequencies 
through synchrotron emission. The same electron pop- 
ulation should also produce hard X-rays via inverse- 
Compton (IC) scattering off cosmic microwave back- 
ground photons. Up to now, however, firm evidence 
for extended non-thermal hard X-ray emission in clus- 
ters was still missing, with a few controversial 11, [l2j 
exceptions, including BeppoSAX [l]| and INTEGRAL 
[li] ] observations of the Coma cluster, BeppoSAX obser- 
vations of Abell 2256 [IH and Chandra observations of 
the Perseus cluster [16]. In addition, under the assump- 
tion of negligible contamination from obscured AGNs, 



[13] reports ~ 2a detections, using the BeppoSAX PDS 
instrument, of hard X-ray non-thermal components from 
a few more merging clusters, including Abell 2142, 2199, 
3376, Virgo and the Ophiuchus cluster itself. 

While ordinary astrophysical mechanisms, including 
merger shocks, can be invoked to explain the non-thermal 
electrons presumably responsible for the observed non- 
thermal activity in galaxy clusters, in the present anal- 
ysis we propose and investigate a novel scenario where 
WIMP annihilations produce, or significantly contribute 
to, said non-thermal population responsible for the hard 
X-ray detection in the Ophiuchus galaxy cluster. In a 
model independent approach, we determine the param- 
eters of the particle DM setups that provide the best 
fits to the INTEGRAL X-ray data, and we compute the 
resulting multi-wavelength spectra. We then compare 
these spectra with the radio data and with the 7-ray 
limits and future prospects for the soon-to-be-launched 
Gamma- Ray Large Area Telescope (GLAST). The high- 
lights of our analysis are: (1) the DM hypothesis will 
conclusively be probed with GLAST; (2) the radio emis- 
sion can in principle also be fitted with the synchrotron 
emission from DM-annihilation- produced e^, as long as 
the average magnetic field in the cluster is of the order of 
0.1 /xG; (3) GLAST might be able detect the monochro- 
matic 7-rays produced in direct DM pair annihilation 
into two photons. 

The flux of e ± produced by WIMP pair annihilations 
depends on the particle DM setup and on the DM den- 
sity distribution. We define a source function Q e (E e ,x), 
which gives the number of e per unit time, energy and 
volume element produced locally in space, as 
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In the equation above, (av)o is the WIMP annihilation 
rate at zero temperature, the sum is over all kincmat- 
ically allowed Standard Model annihilation final states 
/, each with a branching ratio Bf and an distri- 
bution diV/ /dE e , and Apairs^) is the number density 
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of WIMP pairs at a given point x, i.e. the number 
of WIMP particle pairs per volume element squared: 
A/" P airs(£) = Pdm(^)/( 2to dm) ; where p BM stands for the 
DM density. The particle physics framework sets the 
quantity (crv)o, the list of Bf and the mass of the WIMP, 
TO DM . The latter also determines the energy scale of the 
pair annihilation event, and, together with the specific 
final state /, the diV/ /dE e spectral functions, which we 
numerically compute with the Monte Carlo code Pythia 
18]. In addition, Wdjj enters in the determination of the 
local number density of WIMP pairs. 

Once the source function Q e (E e ,x) is determined, the 
spectrum and density are affected by spatial diffusion 
and energy loss processes, usually described - under the 
assumptions of negligible convection and re-acceleration 
effects - by a diffusion-loss equation of the form 
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where dn e /dE e is the number density of electrons per 
unit energy, D is the diffusion coefficient, and 

b(E e , x) = bic + b syn + & C oui + Vcm (3) 

encodes the various energy loss mechanisms [if. 

Knowledge of the distribution of the DM-induced 
population dn e /dE e , of the magnetic field struc- 
ture and strength, as well as of the electron, gas and 
starlight densities allows one to compute the WIMP- 
induced secondary emissions. Specifically, at radio fre- 
quencies the DM-induced emission is dominated by the 
synchrotron radiation of the relativistic secondary elec- 
trons and positrons. IC scattering of the non-thermal 
e ± on target CMB and starlight photons gives rise to a 
spectrum of photons stretching from below the extreme 
ultra-violet up to the soft 7-ray band, peaking in the 
X-ray energy band. Non-thermal bremsstrahlung, i.e. 
the emission of 7-ray photons in the deflection of the 
charged particles by the electrostatic potential of ion- 
ized gas, contributes in the soft 7-ray band. Finally a 
hard 7-ray component arises from prompt emission in 
WIMP pair annihilations, mostly originating from the 
two photon decay of neutral pions, and, at the high en- 
ergy end of the spectrum, from internal bremsstrahlung 
from charged particle final states. The 7-ray spectrum 
extends up to energies equal to the kinematic limit set by 
the WIMP mass, E 1 < m DM , and might feature one or 
more monochromatic lines associated to two-body anni- 
hilation final states where one (or both) of the particles 
is a photon. We refer the reader to Ref. [l| for details on 
the computation of the multi-wavelength emission from 
DM annihilation. 

We show in Fig. [1] the photon flux in the hard X-ray 
and 7-ray bands for three benchmark DM particle model 
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FIG. 1: The hard X-ray and 7-ray spectrum for three D M p ar- 
ticle models, plus a single-temperature MEKAL model [ly] for 
the thermal X-ray emission, compared with the INTEGRAL 
data 0] and with the EGRET upper limit 



accounting for the INTEGRAL data. We also include 
a thermal component from the bremsstrahlung emis- 
sion of the intra-cluster medium, obtained with a single- 
temperature MEKAL model [l9( with the abundance 
fixed to 0.49 compared to the solar value and kT = 8.5 
keV @ . The thermal component was normalized to pro- 
duce the best fit for the INTEGRAL data below 20 keV. 
The DM models were instead normalized to obtain the 
best global fit to the data above 20 keV. We chose DM 
models with Bf = 1 for / = bb, W + W~ , t + t~ , i.e. each 
model pair annihilating into a single Standard Model fi- 
nal state. The three particular final states were selected 
for two reasons: (1) the resulting e spectra dN[ /dE e 
range from the softest (bb) to the hardest (t + t~) possible 
case [l|; (2) the three final states correspond to common 
well-defined cases found in supersymmetric DM models 
2J|. For instance, in the minimal supergravity scenario 
2CJ / ~ bb corresponds to the so-called bulk and funnel 
regions where the neutralino has a relic abundance com- 
patible with the cold DM density, t + t~ is found in the 
coannihilation region (where neutralino pair-annihilation 
proceeds predominantly through scalar tau exchange) 
and W + W~ in the focus point region [21| . This choice of 
benchmark models follows here and generalizes the ap- 
proach of [l| - linear combinations of the considered mod- 
els produce almost any WIMP multi-wavelength emission 
spectrum. 

For each final state, we selected the DM particle mass 
giving the lowest x 2 in the fit to the INTEGRAL data: 
m DM (f = bb) = 46 GeV, m DM (f = W+W~) = 82 
GeV and m,DM(f = t + t~) = 10 GeV. Following Q we 



assumed, for the diffusion coefficient, the form 
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where dg — 20 is the minimum scale of uniformity of 
the magnetic field in kpc and is the average mag- 
netic field in fiG. Notice that, while we take into account 
spatial diffusion in the numerical computation, the effect 
of changing parameters such as Do or ds is minimal on 
the resulting WIMP multi-wavelength spectrum. In the 
computation of the photon flux we neglected the IC from 
starlight and assumed an average thermal gas density 
n = 10~ 3 cm 3 , relevant for the computation of 6coui 
and 6brom in Eq. ©■ In the portion of the spectrum 
shown in Fig. Q] the value of the magnetic field is not cru- 
cial, and was fixed here, for reference, to B = 0.1 fiG. 
As shown in the inset, the fit to the INTEGRAL data 
improves dramatically with the contribution from DM 
annihilation. Also, the emission in the 7-ray band is 
compatible with the EGRET limit [9(, shown for refer- 
ence with a horizontal orange line. 

What is the DM pair annihilation cross section re- 
quired to reproduce the spectra shown in Fig. To 
answer this question we need to integrate the number 
density of DM pairs over the line of sight. In turn, this re- 
quires knowledge of the DM density profile for the Ophi- 
uchus cluster. We follow here the analysis of Ref. [l[, 
and assume, for reference, the DM density profile ob- 
tained in the numerical simulations of Ref. [25| (namely 
p(r) = pog(r/a), with g(x) — exp[— 2(x a — l)/a] and 
a ~ 0.17) and the DM substructure setup outlined in 
[26| . We verified that using other DM profiles changes 
our predictions by less than one order of magnitude. We 
derived from Ref. [53] a virial mass of ~ 1.5 x 10 lb M Q h~ l 
and a virial concentration of ~ 10. We parametrize the 
contribution of substructures via the fraction f s of total 
mass in subhalos, and assume a ratio R s = 5 between 
the concentration parameter in subhalos and that in iso- 
lated halos with equal mass, R s = (c s )/(c v ; r ). In Fig. [2] 
(a) we compute the cross section as a function of the 
DM particle mass, for the three benchmark final states, 
giving the best fit to the INTEGRAL data. We assume 
f s = 0.5 for the lower lines, as suggested by numerical 
simulations [3], while we neglect the contribution from 
substructures for the upper lines. 

While we find rather large values for the pair annihi- 
lation cross section compared to the naive expectation 
(av)o — 3 x 10 -26 cm 3 s _1 motivated by requiring a ther- 
mal WIMP relic abundance compatible with the CDM 
density through simple scaling arguments, the range we 
get is consistent with several examples of supersymmet- 
ric DM models (see e.g. Fig. 15 in []]]). Also, the values 
we obtain for both (av)o and consistent with 

all available particle physics constraints on DM, and are 
compatible with WIMPs being in the right density to- 
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FIG. 2: The preferred pair annihilation cross section (a), the 
integrated 7-ray flux above 100 MeV (b) and the minimal 
branching ratio for the detection of the monochromatic 7- 
ray line at _E 7 = mr>M (c), as a function of the DM particle 
mass. In panel (a) the upper lines refer to the case of no 
substructures, the lower lines refer to the substructure setup 
described in the text, and the gray shaded region is ruled out 
by EGRET [H| and H.E.S.S. [2J] data on the gamma-ray flux 
from the galactic center region. 



day provided, for instance, non-thermal production or a 
modified cosmological expansion rate is assumed at the 
time of WIMP freeze-out [2!|. In addition, uncertainties 
on (1) the dark matter density distribution and (2) the 
galactic and extra-galactic gamma-ray background un- 
dermine the possibility of ruling out WIMP models with 
large pair-annihilation cross sections via gam ma-ray data 



from the center of the Galaxy ,2 
galaxies [3 



2J, l3fj, l3l|^iearby 
33j, |34 j , and from the galactic halo [35| . For 



instance, in the upper panel of fig. [2] we shade in gray 
the region ruled out by EGRET and H.E.S.S. data from 
the galactic center region, assuming a cored dark mat- 
ter profile following the analysis of Ref. 3(| 31 1 . Clearly, 
gamma-ray data from the galactic center do not rule out 
the range of cross sections we find. 

Similarly, antimatter constraints not only depend on 
the dark matter density across the hole galactic halo, but 
are also affected by sizable uncertainties in the diffusion 
and propagation of cosmic rays in the Galaxy [35[ . Even 
for cross sections as large as to account for the EGRET 
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FIG. 3: The spectral energy distribution for the multi- 
wavelength emission of three DM particle models and for the 
thermal X-ray emission. We also show the INTEGRAL 
and radio data. 



data on the galactic gamma-ray emission [36[ (hence at 
the level of the cross sections in the shaded gray region in 
fig. [3]), antimatter fluxes from WIMP annihilations are in 
general compatible with available data [37| • Finally, both 
the search for energetic neutrinos from WIMP annihila- 
tion in the Sun or the Earth, and direct WIMP detection, 
depend on the WIMP-nucleon scattering cross section, a 
quantity which is unrelated to the WIMP pair annihila- 
tion. Even in special models, such as supersymmetry, the 
range of predictions is so wide to make it impossible to 
constrain (av)o with neutrino fluxes or direct detection 
searches 38|, [39j. In conclusion, the WIMP pair anni- 
hilation rates we consider here, while larger than what 
naively expected, are compatible with WIMP cosmolog- 
ical DM production and DM searches. 

Panel (b) in Fig. [5] shows the integrated 7-ray flux 
above 0.1 GeV for the best fit models as a function of 
m DM ■ In all cases we find that the expected 7-ray flux is 
well above the anticipated GLAST LAT integral flux sen- 
sitivity, estimated to be around a few xl0~ 10 cm~ 2 s _1 
[40l | . Fig. [3(c) shows the branching ratio (crv) to t / (<^ v ) 77 
for the monochromatic DM DM — ► 77 channel needed to 
obtain, for the best fit models, the detection of at least 
10 photons 1 with E 1 — m DM . Notice that the values 
shown are independent of the assumed DM profile. While 
(c^)tot/ (°"i>) 77 is entirely model dependent, the range we 
obtain is generically consistent with what is expected e.g. 
in supersymmetry [4 11 ] , and especially in next-to-minimal 




1 The branching ratio corresponding to a larger photon flux can 
be obtained by linearly rescaling the lines in Fig. l2l (c). 



FIG. 4: Preferred values of the cluster magnetic field as a 
function of the DM particle mass, for three representative 
pair annihilation final states. 



supersymmetric extensions of the Standard Model [42j. 
GLAST can therefore easily detect a sizable number of 
monochromatic energetic 7 rays, depending on the spe- 
cific DM particle model. 

We illustrate the whole broad-band spectrum for the 
DM annihilation interpretation of the INTEGRAL hard 
X-ray emission in Fig. [3] Particularly crucial here is the 
question of whether the radio data on the Ophiuchus clus- 
ter are compatible with our predictions. This question 
depends on the assumed value of the average cluster mag- 
netic field B. In the Figure we use values of B giving the 
best fit to the radio data, namely B(f = bb) = 0.15 /iG, 
B(f = W+W~) = 0.1 plG and B(f = t+t~) = 0.18 fiG. 
Clearly, a linear superposition of the various final states 
can yield a remarkably good fit to the available radio 
data. We also indicate the anticipated GLAST sensitiv- 
ity, which shows that a sizable flux of 7-rays is expected 
if the DM annihilation scenario is indeed correct. 

As pointed out in the analysis of Ref. [6( , simple IC in- 
terpretations of the INTEGRAL data force the estimate 
for the magnetic field in the range B ~ 0.1 — 0.2 /iG, 
which agrees with the present analysis. While this range 
is compatible with the value obtained from the hard X- 
ray emission in Coma [l3j ]. it is below the estimates ob- 
tained, for other clusters, through Faraday rotation mea- 
sures [43| (see however the discussion in (lal). We show 
in Fig. 2] with black and red lines the values of B respec- 
tively giving the best fit to radio data and exceeding by 
5-cr the measured radio flux in at least one bin. While 
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values of B above the red lines are disfavored by radio 
data, they are not strictly ruled out if one accounts, e.g., 
for a radial dependence of the magnetic field [l[ . 

The detection of 7-rays from cosmic rays in clusters 
of galaxies is g enerically expected to be possible with 
GLAST 0,0- Specifically, the analysis of Ref. in- 
dicates that the Ophiuchus cluster will be detectable by 
GLAST provided the ratio of the energy density of cos- 
mic rays to that of the thermal gas is larger than 0.3% 
to 8.8%, depending on assumptions on the cosmic ray 
spectral index and radial distribution. This warrants a 
systematic comparison between the expected 7-ray flux 
from cosmic rays and from DM in clusters, which is cur- 
rently under way [4(1 ]. In general, though, the 7-ray 
spectra expected from cosmic rays (a power law above 
~ 1 GeV, with a spectral index depending on the pri- 
mary cosmic-ray protons spectral index) are significantly 
different from what expected from DM annihilation, see 
Fig. [T] and [3] If DM annihilation is responsible for most 
of the hard X-ray emission in the Ophiuchus cluster, 
GLAST will collect such large statistics that discrimina- 
tion from a cosmic ray 7-ray emission appears reasonably 
feasible. 

In summary, we showed that the origin of the non- 
thermal particles plausibly responsibly for the recently 
firmly discovered non-thermal hard X-ray emission from 
the Ophiuchus cluster might be generated by electrons 
and positrons produced in WIMP pair annihilations, pro- 
vided the rate for the latter is large enough. This sce- 
nario is compatible with all observational information on 
the cluster, with particle DM production and searches, 
and, more importantly, will be thoroughly tested by 
GLAST. The future 7-ray telescope might even detect 
the monochromatic two-photon emission provided the 
particle model has a large enough branching ratio in that 
channel. Other galaxy clusters exhibiting non-thermal 



activity [17| will also be outstanding sites for GLAST to 
look for signatures of WIMP DM pair annihilation [4(| . 
Finally, radio data from the Ophiuchus galaxy cluster can 
also be accounted for in the DM annihilation scenario, as 
long as the average magnetic field in the Ophiuchus clus- 
ter is below « 0.2 /xG. 

The author wishes to thank Tesla Jeltema for several 
insightful comments and for help with the computation 
of the thermal X-ray emission. 
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